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ON  SIMPLE  FINITE  CAPACITY  DAMS  AND  QUEUES 


1.  Introduction 

V*  are  concerned  with  various  properties  of  the  following  two 
■ode la,  which  arise  in  the  theory  of  daas  and  of  queues.  In  the  first 
example  we  consider  a  daa  of  finite  capacity  b  (a  positive  real 
nuaiber)  fed  by  inputs,  where  the  else  of  each  being  a  random  variable 
having  a  negative  exponential  distribution  with  mean  pf1  ,  which  occur 
In  a  time-homogeneous  Poisson  process  with  parameter  A(0  <><■>)  . 

The  release  of  water  is  continuous  at  unit  rate  per  unit  time  unless 
the  dam  is  empty,  and  any  part  of  an  input  which  raises  the  total  content 
b  overflows  and  is  lost. 

In  the  second  example  we  assuae  customers  arrive  at  a  counter 
with  a  single  server  in  a  time -homogenous  Poisson  process  with  parameter 
A  (0  <  >  <  •)  .  Customers,  who  are  served  according  to  a  first-come 
first-served  queue  discipline,  have  independently  and  identically  dis¬ 
tributed  negative  exponential  service  tiats  with  mean  p.  1  .  There 
is  deterministic  customer  impatience,  so  that  if  a  customer  arrives  to 
find  that  he  has  a  wait  exceeding  b  before  coomnclng  service,  then 
he  departs  forever  from  the  system  without  being  served. 

Each  of  the  above  two  models  have  been  discussed  by  a  number  of 
writers.  Phatarfod^  has  made  use  of  Wald's  identity  to  obtain  the 
distribution  of  the  time  to  *lrst  emptiness  (with  and  without  overflow 
being  allowed)  in  the  daa  model;  Veesakul  and  Yeo  ^  have  obtained 


these  and  the  time -dependent  dam  content  distribution.  In  the  case 


of  the  queue  various  similar  results  for  waiting  time  and  queue  size 


have  been  obtained  by  Barrer 


[1] 


^  and  Tzkdcs^  and  othere. 


Finch 
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It  la  the  purpose  of  this  note  to  show  how  many  of  these  quantities 
can  he  obtained  directly  from  the  appropriate  Kolmogorov  equation  by 
forming  and  solving  second  order  linear  differential  equations  with 
constant  coefficients.  This  gives  a  uniform  method  for  finding  these 
and  other  similar  quantities. 

We  first  discuss  the  problems  of  the  dam  model  and  then  the 
queueing  model;  the  former  is  rather  simpler  to  deal  with,  as  there 
cure  fewer  complications  beyond  the  barrier  b  than  exist  in  the 
latter  case. 

2.  Dams 

We  denote  the  dam  content  at  time  t  by  Z(t);  we  are  Interested 
in  the  behavior  of  Z(t)  as  time  passes.  We  discuss  several  quantities 
individually . 

2.1  First  Baptlness  Before  Overflow 

We  define  the  probability  of  first  emptiness  by  time  t  before  over¬ 
flow  given  the  initial  content  by 

G(t  ,  c)  -  Pr  (Z(t)  -  0  for  scat  t  £  t  ; 

0  <  Z(v)  <  b  ,  0  <  v  <  t  |  Z(0)  -  c  ,  0  <  c  <  b) 

and 

00 

0(©  ,  c)  -  J  e_atG(t  ,  c)dt  R1  ©  >  0  . 

t-0 

Considering  the  time  interval  [0  ,  t  ♦  6)  for  6  >  0  in  two 
phases  [0  ,  6)  and  [6  ,  t  ♦  6)  we  have 
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G(t  ♦  6  ,  c)  -  (1  -  >fc)G(t  ,  c  -  6) 


♦  Xt  J  G(t  ,  u  ♦  6Hi«"U^U'c)du  ♦  0(6)  0  <  t  <  oc  , 

u«c 

which  yield*  the  backward  Kolmogorov  equation 

-  X0(t  ,  c)  ♦  V J  G(t  ,  u)e"^(c"u)du  .  (2.1) 

UvC 

TWking  train  forma  with  respect  to  time  we  find  that 

b 

0* (0  ,  c)  ♦  (*  ♦  0)0(0  ,  c)  -  >ii e^CJ  0(0  ,  u)e"*iUdu  (2.2) 

u-c 

where  differentiation  is  with  respect  to  c  .  Ve  now  differentiate  (2.2) 
with  respect  to  c  and  obtain 

0"(O  ,  c)  ♦  (>  ♦  0  -  u)0'(O  ,  c)  -  u©0(©  ,  c)  -  0  ,  (2.3) 


which  is  a  second  order  linear  differential  equation  with  constant  coef¬ 
ficients.  The  boundary  conditions  are  found  by  inspection  of  the  process 
to  be 


0(0  ,  0)  -  0'1 

0 ' (0  ,  0)  -  lia  0 ' (0  ,  c)  -  -  (>  ♦  0)0(0  ,  0) 
o*0 

b 

♦  >iij  0(0  ,  u)e"uudu  . 

u»0 

The  solution  of  (2.3)  1*  readily  obtained  by  standard  methods 
and  is 

£1C  £2C 

0(0  ,  c)  *  a^e  ♦  o^e  0  <  c  <  b 


where  may  be  functions  of  0  but  not  of  c  ,  and 


(2.4) 


(2.5) 
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(2.6) 


*1,2  "  -  ♦  »  -  u)  ±  4((x  ♦  ©  -  »Y ♦  4%}* 

.  -  $(X  ♦  0  -  u)  ±  £((X  ♦  ©  ♦  ii)2  -  4X*i)^ 

Use  of  the  boundary  conditions  (2.4)  enables  us  to  determine  and 

Og  ;  after  some  straightforward  algebraic  manipulations  we  obtain 

(So-U^iC  (t1-|i)b'*C9c 

(C,  -  U)e  ^  -  (Co  -  u)« 

1,(8  • c)  *  - tr^ - r?r^r7B —  (2-7) 

#((«!  -  u)«  “  -  (t2  -  ll)«  1  ) 

which  agrees  with  (J.l)  of  [5]  and  with  (2.6)  of  [3].  The  result  (2.7) 
may  be  inverted  to  give  G(t  ,  c)  ;  this  has  been  done  in  [3],  and  the 

result  Is  in  terms  of  Bessel  functions.  Similarly  all  other  results  in 

this  note  for  time  transforms  of  content  and  waiting  time  may  be  inverted. 

Ve  note  that  the  probability  G(»  ,  c)  that  emptiness  occurs  before 
overflew  is  given  by 

-^(b-c)-Xc  -Xb 


G(.  ,  c)  -  lie  00(0  ,  c)  - 
0-0 


Xe 


-Vtb 


-  tie 


2.2  First  Emptiness  When  Overflow  JSML  Occur 

In  the  previous  problem  the  process  ended  when  either  the  dam  emptied 
or  overflowed;  we  now  continue  until  emptiness  occurs,  regardless  of 
whether  there  has  been  overflow  or  not.  Only  a  slight  modification  of 
the  boundary  conditions  is  required.  We  define 

G  (t  ,  c)  «  Pr  (Z(t)  -  0  for  seme  T  £  t  ; 

Z(v )  >  0  ,  0  <  v  <  t  |  Z(0)  -  c  ,  0  <  c  <  b) 


A©  ,  c)  -  J  e'®tG#(t  ,  c)dt 
t-0 


R1  0  >  0  . 
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The  backward  Kolmogorov  equations  and  the  differential  equations 
are  of  an  identical  fora  to  (2.1)  -  (2.)),  so  that  the  solution  is  of 
the  fora  (2.5 );  however,  the  boundary  conditions  are  now 

0*(O  ,  0)  -  O'1 

0*(O  ,  0)  -  lUn  0*'(©  ,  c)  -  -(X  ♦  0)0*0  ,  0) 
c-0 

♦  >4 J  0  (0  ,  uje'^du  +  X0  (0  ,  bje”^3 

u*0 


Proceeding  as  above  we  obtain 


(^.-uJb^c 

#  (ex  -  n  ♦  X)e  ^  -  (*2  -  ^  ♦  A)e 

0  (» , c) - —  ' 

«((*2  -  H  ♦  *)e  -  -  U  ♦  X)e  1  } 

which  agrees  with  (5.1)  of  [5]  and  with  (2.10)  of  [5].  We  note  that 

lia  0'V(0  ,  c)  -  1,  so  that  eventual  emptiness  is  certain. 

©~0 


2.3  The  Duration  of  a  Wet  Period 

A  wet  period  is  the  time  from  the  entrance  of  am  input  into  an 
empty  dam  to  the  next  moment  the  dam  again  empties;  we  define 


G(t)  -  Pt(Z(t)  *  0  for  same  t  £  t  ; 

Z(v)  >0  ,  0  <  v  <  1  J  input  at  time  sero) 

0(0)  -  J  e'®tG(t)dt  R1  0  >  0  . 

t-0 

Obviously, 


(2.8) 


(2.9) 
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♦  M-' 


l  c 

0(0)  -  J  ,  c)dc  ♦  0*(O  ,  b)i 


.1* 


c-0 


(k-u)* 

($1  -  nKfcg  -  |x  ♦  X)e  -  •  4)(ftx  -  4  ♦  x)e  ^ 

T?]^nr 


(2.10) 


«*(*!  -  ^(Cg  -  4  ♦  X)e 


TttFT 
-  (Cx  -  n  ♦  ^  ) 


2.4  The  Content 

We  define  the  daa  content  distribution  at  tlan  t  given  the 
initial  content  by 


’(x  ,  t)  -  Pr(Z(t)  £  x  |  Z(0)  -  c  ,  0  £  c  £  b) 


and 


*(x  ,  0)  -  J  e"®V(x  ,  t)dt  Rl  0  >  0  . 

t-0 

In  this  case  ve  sake  use  of  the  fonmrd,  rather  than  the  backward, 
Kolaogorov  equation.  Ve  have  [5 ] 

§jF(x  ,  t)  -  |*(x  ,  t)  -  -  >T(x  ,  t)  >.»' 

♦  V  J  F(u  ,  t)e~^x”u^du  0  £  x  <  b  (2.11) 

MmO 

F(x  ,  t)  -  1  x  ^  b 


and  consequently 

t'(x  ,©)-(>•*■  ©)*(x  ,  0) 

-  -  Hie”^*  J  f(n  ,  Oje^du  -  U(x  -  c)  , 

u»0 

where  U(x)  is  the  unit  step  function.  Differentiating  with  respect 
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to  x  ve  obtain 


f"(x  ,  0)  -  (X  +  0  -  *  ©)  -  U©f(*  »  ©)  s  *  ^U(x  -  c) 

0<x<b,x*c, 


(2.12) 


with  the  boundary  conditions 


f(b  ,  0)  -  0_1  ;  f(0  ,  0)  «  11a  f'(x  ,  0)  =  (X  +  0)*(0  ,  0) 

x-*0 

♦(c*  ,  0)  -  *(c-  ,  0)  ;  *'(c+  ,  0)  -  *’(c-  ,  0)  -  1  . 


(2.13) 


New  aquation  (2.12)  baa  the  solution 

*2* 


*(x  ,  0) 


r  v  x  vgx 

|  V  +V 

<|  V  X  V  X 

[v  ♦  V  ♦« 


0  Si  *  £  c 

C  ^  X  ^  b 


(2.14) 


where  ai  »  ’  ^1  ’  ^2  do  not  depend  on  x  and 

vi  2  =  £(*  ♦  ®  "  n)  *£((*♦©♦  M-)2  -  4>p.}^  • 


(2.13) 


We  determine  a1  ,  a?  ,  b^  ,  bp  from  (2,13)  and  (2*llf)»  aft«r  some 
algebraic  manipulations  ve  obtain  (c.f.  (6.1)  of  [3]) 
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-(vl+v2)c  ^  -vl(b-c)  ~v0(b-c ) 


Jv1-0)e 


-V  (b-c)  'l 
■(Vg-O)e  2  | 


{(^.v2,.V(^.Vi)eV} 


[  -v  b  -v  b 

j/(x,Q)  =  <,  x  ^©(v^VgilO+S-Vgk  2  -{>+«-v1)e  |  J 


-1 


0  1  x  £  c 


f  'V  -V 1 

|(v2-0)(>4«-v1)e  x  +(Vl-8)(^-v2)e  *  j 


-v  (b-x)  -v  (b-x) 
x  le  -e 


{< 


-v  b 

xi(^-v2)e  -(X+3 


}[#(vV 

-V^b}j 


-1 


c  £  x  £  b 


The  at*  denary  dam  content  distribution  nay  be  obtained  directly 
from  the  forvard  Kolmogorov  equation  or  aa  the  limi  t  as  t  -•  <*  of  the 
time -dependent  distribution.  If  F(x)  *  ^lm  F(x  ,  t)  then 

-(n-X)x 


F(x)  =  11m  0 fix  ,  0) 
O-iO 


U  -  Xe 
p.  -  Xe 


r 

1 . 


-  Cm-"X  )b 


°£*£b 

x^b 


3.  Queues 

3.1  First  Idleness  of  the  Server 

If  we  consider  first  idleness  of  the  server  without  the  virtual 
waiting  time  exceeding  b  then  the  problem  is  identical  to  that  of 
first  emptiness  without  overflow  in  the  dam  model  We  thus  restrict 
ourselves  to  the  case  where  the  waiting  time  nay  exceed  b  ,  but  no 
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customer  will  accept  a  wait  greater  tnan  b  . 

We  let  Y(t)  be  tbe  virtual  waiting  tine  at  time  t  .  We  define 

H(t  ,  c)  >  Pr  (Y(t)  *  0  for  seme  t  £  t  ; 

Y(v)  >  0  ,  0  <  v  <  t  |W(0)  =  c) 

/  00  -®t 

•(0  ,  c)  -  J  e  H(t  ,  c)dt  R1  0  >  0  . 

t-0 

The  backward  Kolmogorov  equation  for  this  quantity  is 


-  XH(t  ,  c) 


^H(t  ,  c)  +  ^fl(t  *  c) 


<  ♦  Hi  j  G(t  ,  u)e‘U^U'C^du  0  <  c  <  b 


UsC 


c  >  b, 


so  that  by  taking  time  transforms 


>♦(0  ,  c) 


CO 

»’(©  ,  c)  +  9«(9  ,  c)  =  <  +  Hie4Cj  «(9  ,  uje^du  0  <  t  <  b 


Use 


c  >  b. 


Hence 
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©"(©  ,  c)  ♦  (X  ♦  ©  -  n)*'(©  ,  c)  -  n©*(©  ,  c)  o  0  0  <  c  <  b 

(3-D 

®'(©  ,  c)  +©*’(©,  c)  •=  0  c  >  b  , 


with  the  boundary  conditions 

•(©  ,  0)  -  O-1  ;  ♦(©  ,  b+)  -  ®(©  ,  b-) 

♦  ’(©  ,  0)  -  11a  ♦*(©  ,  c)  ■  •  (>  ♦  ©)*(©  ,  0) 

o*0  (3.2) 

♦  \i J  •(©  ,  uje^du  . 

u«0 


This  is  very  similar  to  the  case  of  first  eaptiness  in  a  daa  with 

an  extension  to  c  >  b  ;  we  obtain  in  this  case 

(l«-u)b4€  c  (i-iOb^c 

Ui-uKt,**)*  1  -(^-u)(^)«  1  ^ 

•(0»c) - - rr^ —  (5*3) 

0((41-4)(42*Hi)e  ^  -(^-U)^)*  1  ) 

where  ^  ^  are  given  by  (2.6)  . 


3.2  Iteration  of  Busy  Period 

A  busy  period  lasts  fraa  the  point  of  arrival  of  a  customer  at  a 
free  counter  to  the  next  instant  a  departing  customer  leaves  the  server 
idle;  we  define 


H(t)  «  Pr  {  Y(i)  «=  0  for  bom  t  £  t  ; 

Y(v)  ^  0  t  0  <  v  <  t  |  a  custoaer  arrives  at  tiae  aero) 


♦(©)  -  J  e’®tH(t)dt  R1  0  >  0  . 

t-0 

It  follows  immediately  froa  the  results  of  (3.1)  that 
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(3-  M 


w 

•(©)  -  J  4a  "^©(O  ,  c)dc 


c*0 


(i,  -  U2  " 

\l(i±  *  ©)e  -  4(^2  +  0)e 

crrnrnr 

-  (t,  -  u)(^  +  «)•  1  } 


d2  -  u)f 


oux  -  n)(e2  ♦  «)• 


3.3  The  Waiting  Time 

We  define  the  time-dependent  waiting  time  distribution  as 
V(x  ,  t)  -  Pr  (Y(t)  £  x  |  Y(0)  =  c  ,  0  £  c  £  b} 

and 

ft(x  ,  0)  *  J  e‘®tW(x  ,  t)dt  R1  ©  >  0  . 

t-0 

Proceeding  as  for  the  dam  content  distribution  ve  obtain  the 
following  equations: 


-  >W(x  ,  t) 

|^W(x  ,  t)  -  |©(x  ,  t)  «  <♦  >4 J  *W(u  ,  t)e“u(x'u)du  0  <  x  <  b 

UaO 

-  Ml(b  , 

*  >41 J  W(u  ,  t)e"*^x_uW  x  >  b 

u*0 
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(3.5) 


0"(x,©)  -  (Vtf-uXl'U,®)  -  H«l(x,0)  -  -^U(x-c)  0<x<b  ,  x  {  c 

n"(x,o)  +  (p.-oXi'U,©)  -  ^©o(x,o)  =  -4  x>b  , 

vith  the  boundary  conditions 

n(c+,  0)  *  n(c-,  0)  ;  n'(c+,  ©)  «  a'(c-,  e)  -  1 

0(b+,  0)  =  fl(b-,  0)  ;  ft’(b+,  ©)  «  £l'(b-,  ©) 

0(-  ,  «)  *  o'1  ;  fl’(0,  ©)-(>♦  0)0(0  ,  ©)  . 

The  solution  to  (3*5)  1* 

0  £  *  £  c 

c  £  X  £  b  , 

^  r.e"4*  ♦  O-1  x  £  b 

where  »  c2  *  \  40  not  <**pend  on  x  and  g  Is 

given  by  (2.15).  Together  vith  (3*6)  ve  obtain  frca  (3-7)  that 


V  X  V  X 

1  2 

Cl*  +  C2* 


i  V  V2X  -1 
fi(x  ,  ©)  -  <,  V  ♦  V  ♦  ® 


(3-6) 


(3-7) 
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n(* 


(Vv2)c 


f  -v  (b-c)  -v  (b-c)  l 

|(v1-0)(n+v2)e  -(Vg-oKu+v^e  j 


V 


*  (W-v9)e  A  -(>4®. 


,-l 


X  [^(v^Vg)  |(v1-Hi)(X4«-Vg)*  2  -(v2-»u)(>4«-v1)e  1  |J  °  £  x  <; 


if  -v.c  -v  c  'i 

,0)  -  <  |-(v2-0)v>4«-v1)e  1  +(Vl-0)(>4«-v2)e  2  | 


f  -v  (b-x)  -v  (b-x) 'j 

|(v1-*u)«  -(Vg^u)®  I 


[o(Vv2)  {(V1-Hi)(>v4«-V2)e  2  -(Vg^iOO^-v^e  2  }] 
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♦  o 


0  £  x  £  b 


We  note  that  W(0  ,  t)  and  1  -  W(b  ,  t)  are  respectively  the 
probability  of  a  customer,  vho  arrives  at  time  t  ,  commences  service 
insedlately  and  is  never  sei-ved;  ve  have  the  time  transforms  of  these 
quantities. 

As  for  the  dam  model  the  stationary  waiting  time  may  be  obtained 
directly  from  the  fonrard  Kolmogorov  equation,  or  from  the  time -dependent 
solution  as 


W(x)  «=  lim  W(x  ,  t) 

t-«oo 


TlTTJb 


0  £  x  £  b  . 
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